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We investigate coherent and incoherent tunnehng phenomena in conditions of crossing diabatic 
potentials. We consider a model of two crossing parabolic diabatic potentials (left (L) and right (R)) 
with an independent of coordinates constant adiabatic coupling Ui2- As a result of the coupling and 
level crossing avoiding, we get the asymmetric double-well lower adiabatic potential with a variable 
shape depending on a value of a continuous parameter b (which describes in the limit 6=1 two 
identical parabolic diabatic potential crossing and in the limit fe — > 00 one-well and linear diabatic 
potentials crossing) . We show that the doublet structure of levels (generic for double- well potentials) 
is remained valid as long as the transition matrix element Hlr (or tunneling splitting) is smaller 
than characteristic inter-level spacings Ar (the latter ones in own turn decrease upon increasing of 
a difference between the diabatic potential minima 5Elr)- We calculate the non-adiabatic factor, 
i.e. Hlr as a function of Ui2- In the diabatic limit {U12 —> 0) Hlr goes to zero, and in the adiabatic 
limit {U12 — > 00) the tunneling transitions do not depend on the upper potential. In the over-barrier 
energy region Hlr is an oscillating function of U12, due to the resonances between the states in the 
lower and in the upper adiabatic potentials. In the case Hlr > Ar, any level from the shallow 
L-well is coupled by the tunneling to several levels in the ii-well, and the transitions lose their 
coherence. A new phenomenom emanated from this oscillating dependence of Hlr on U12, namely, 
multiple coherent - incoherent regime transitions for the upper adiabatic potential state evolution, is 
our main concern in this paper. The problem is not only of intellectual interest but also of relevance 
to various molecular systems undergoing conversion of electronic states or isomerization reactions. 
Our model exhausts all cases practically relevant for spectroscopy of non-rigid molecules, and can 
capture many of the features exhibited by experiment. 

PACS numbers: 05.45.-a, 72.10.-d 



I. INTRODUCTION 

Double-well potentials appear in various contexts in physics and chemistry. For example the simplest pattern 
of almost any molecular reactive system (with two stable configurations identified as a reactant and as a product) 
corresponds to the model potential energy formed by two multidimensional more or less parabolic surfaces shifted 
relative to each other. Although the ID asymmetric double-well model is idealized, it can be very useful for a 
qualitative discussion to gain more insight into complex multidimensional dynamic molecular properties for which 
exact or even approximate theoretical results are not available, thus throughout what follows we will consider ID case 
only. 

In the classical limit for the energy region E < Ub (where Ub is the potential barrier, separating the both, say L and 
R wells) which will be referred further on as the tunneling region, the behavior in the both wells are fully decoupled 
and therefore one well is independent from other. As it is a common wisdom nowadays in quantum mechanics even 
for E < Ub the particle can tunnel between the wells. It admixes the L and R well localized states, thus allowing 
an under-barrier tunneling mechanism. The extent of this delocalization is larger in the states close to the top 
of the barrier, and it is maximal when the unperturbed levels on the opposite sides of the barrier are degenerate 
(the reason is immediately clear by looking at the standard textbook expressions for the tunneling probability and 
splitting [1]). For the symmetric case this tunneling level splitting leads to coherent quantum oscillations typical for 
any two-level system. For asymmetric double well potentials, pairs are not in coincidence any more, the tunneling is 
suppressed, except for certain critical values of model parameters for which the levels are brought in resonance again, 
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and the problem becomes more tricky. We have recently shown [2] that one can successfully attack this problem 
by a semiclassical solution of the Schrodinger equation for ID asymmetric double well potential with one-parameter 
dependent shape 



Ui is written in dimcnsionlcss form measuring energy in the characteristic frequency of the oscillations say around 
the left (L) minimum Qq, and the coordinate X is measured in the units ao of the inter- well distance (and in "God 
given unit" we put h = 1, except where explicitely stated to the contrary when it is necessary for understanding or 
estimations). The climensionless parameter b allows us to change the shape of the R well and to consider both limiting 
cases, namely a traditional symmetric double-well potential (6 = 1), and a decay potential for 6 ^ oo. Behavior in 
the latter limiting case is also well known, there is a continuum spectrum of eigenstates for X — > oo and incoherent 
decay of quasi-stationary states from the L - well. 

For a small asymmetry 6 ~ 1 the ground state doublet Eq in the potential (1.1) 



where Hj^^ is the ground state splitting for the symmetrical double well potential. 

It was shown in [2] that for any asymmetric double-well potential the behavior depends crucially on a dimensionless 
parameter A that is, roughly speaking, a ratio of characteristic frequencies for low-energy in-well oscillations and 
inter- well tunneling. For A ^ 1, there are well defined resonance pairs of levels, and so-called survival probability 
(i.e. the probability for a particle initially localized in one well to remain there) has coherent oscillations related to 
resonance splitting. However, for A — > oo for any finite time scale, there are no oscillations for the survival probability, 
and there is almost an exponential decay with the characteristic relaxation time cx determined by Fermi golden 
rule. In the NMR language this relaxation time can be associated to the so-called dephasing T2 time. Thus one can say 
that tunneling destroys coherent behavior and it can be associated with dephasing processes in the phenomenological 
Bloch theory of quantum relaxation. Explicitely for asymmetric double-well potentials 



where An is a typical level spacing for the final states. In the case A 1 one may not restrict himself to the only 
resonance pair levels. The number of levels perturbed by tunneling grows proportionally to \/A, or, by other words, 
instead of isolated pairs there appear the resonance regions containing the sets of strongly coupled levels. At the 
intermediate values of A ~ 1 one has a crossover between both limiting cases, namely the exponential decay with 
subsequent long period recurrent behavior (more longer the larger A). In this region, A ~ 1, and in the vicinity of 
quasi-stationary states of the L- well the doublet-like eigen-spectrum turns into the equidistant one. 

It is particularly instructive to look to this result from a slightly different point of view related to a striking and 
still enigmatic phenomenom - quantum chaos. Perhaps the first successful quantitative criteria relating the classical 
ergodic theory to quantum molecular dynamics was formulated long ago by von Neumann and Wigner [3]. According 
to [3] a system has the ergodic behavior if it has: 

(i) equidistant spectral distribution (i.e. no degenerate states); 

(ii) time decay of correlations for any observable. 

Let us emphasize that the criteria claim that quantum manifestations of the classical chaos are related to specific 
spectral features of a system under consideration, and not to some kind of its irregular time evolution. From the first 
sight it might seem that ergodic behavior is impossible in one dimensional systems which are integrable in classical 
mechanics. However, in fact our model examples should be considered as two classically decoupled systems {L and R 
wells) interacting only via quantum tunneling. One can call this phenomenom as tunneling induced ergodicity, and 
for more details, see e.g. [2], (and also, [4], [5]). It turns out that the time evolution of the localized initial state is 
governed by the interplay of two physical parameters, recovering period time and decay rate. The both parameters 
depend on the spectral representation (or distribution function) of the initial state. It is shown in [2] that the both 
listed above criteria are satisfied for strongly asymmetric double-well potentials where highly excited states in the R 
well are strongly perturbed by tunneling from the shallow L well. The condition to have this phenomenon (tunneling 
induced ergodicity of final states) reads in our notation as 




(1.1) 




(1.2) 



\Hlr\>Ar. 



(1.4) 
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In this paper wc will show that similar phenomena can take place also for dclocalized initial states, where non-adiabatic 
coupling gives rise chaotic behavior. It will be referred in what follows as non-adiabatic transitions induced ergodicity. 

Of course an isolated double-well potential is only an idealization of any real molecular system. The applicability 
of such an idealization must be decided separately for each system or process in question. However even in the cases 
where such a model is not justifiable the calculations we performed are nonetheless instructive. But in this paper 
we make one step further. In a typical problem of chemical dynamics or molecular spectroscopy, the double-well 
potentials can appear as a result of level crossing phenomena, and the consideration of only the isolated double-well 
potential (lower adiabatic potential) can be justified only if the gap occurring in the spectrum at the avoided level 
crossing point is much larger than all other characteristic energy scales of the problem. However, evidently it is not 
the case for example if we are interested in the calculation of vibrational - tunneling spectra of non-rigid molecules, 
or reactive complexes with more than one stable configuration. The lowest multi-well potential of such systems is 
formed from one well diabatic potentials crossing corresponding to each stable configuration. Apart from the lowest 
potential, the upper adiabatic potential with its minimum above the maximum of the lowest potential should be also 
taken into account for these situations (see Fig. 1). In the most of the calculations of tunneling splittings in the 
ground and low excited vibrational states the coupling to the upper potential are neglected, what is certainly correct 
only for strong enough adiabatic coupling. The same situation takes place for systems undergoing the Jahn - Teller 
effect, where the interference of the diabatic states occurs [6]. In all these situations the adiabatic coupling removes 
diabatic level crossing, and the diabatic levels are replaced by the adiabatic ones. Let us repeat that only in the case 
of a large adiabatic splitting one can restrict oneself to the only lower adiabatic potential and neglect any influence of 
the upper adiabatic potential. However, in a general case of arbitrary adiabatic splittings, intra-well and inter-wells 
dynamics depends on the both adiabatic potentials (i.e. on tunneling splittings and adiabatic interactions). 

In the fundamental problems of chemical dynamics and molecular spectroscopy, the transitions from the initial to 
final states can be treated as a certain motion along the potential energy surfaces of the system under consideration. 
These surfaces in own turn are usually determined within the Born - Oppenheimer approximation. However, the 
approximation becomes inadequate for the excited vibrational states, when their energies are of the order of electronic 
inter level energy spacing or near the dissociation limit. In the both cases the non-adiabatic transitions should be 
taken into account, and the most of the non-radiative processes occur owing to this non-adiabaticity. The typical 
examples investigated in the monography [7], are so-called pre-dissociation, singlet-triplet or singlet-singlet conversion, 
and vibrational relaxation phenomena. 

To treat this kind of level crossing (Landau - Zener (LZ)) problems usual textbook consideration utilizes the outset 
within a limited electronic subspace which is completely spanned by a finite set of Born - Oppenheimer or adiabatic 
electronic states. However, because these states obey the noncrossing rule it may be desirable technically to transform 
the states into the diabatic representation in which the diagonal matrix elements of the electronic Hamiltonian in the 
subspace can cross, and the off-diagonal interactions appear as scalar coupling potentials. 

The major concern of this paper is with the construction and solution of a model for two asymmetric diabatic level 
crossing phenomena. The rest of our paper begins in section II with a formulation of our model and with a discussion 
of basic methodical details necessary for our study. Section HI contains our main results. We derive the criteria for 
reversibility and coherent or incoherent tunneling for crossing diabatic potentials. Our conclusion section IV deals 
with miscellaneous subjects related to the diabatic level crossing phenomena. 



II. MODEL POTENTIAL AND BASIC RELATIONS 



As a model for diabatic potentials in this paper we choose two non-equivalent parabola 

C/i = i (1 + X)^ C/h = i (l - 2X -f ^) (2.1) 

with a symmetrical crossing in the point X = 0. Upon increasing the well asymmetry 

SELn = -^ (2.2) 

the potential Ub. is converted from a simple parabola at 6 = 1 to a linear potential at 6 ^ oo. Owing to the adiabatic 
coupling Ui2 (which we assume for simplicity independent of coordinates) we get the lower double- well and the upper 
one- well adiabatic potentials (see Fig. 1). 

At arbitrary values of the parameters U12 and b to find eigenstates and eigenfunctions for our model potential we 
should solve the coupled Schrodinger equations 

- ^ + iHUl{X) - E)@L = 1^U^2Qr ; + i\Ur{X) - E)Qn = i'U,2@l , (2.3) 



4 



which can be written as one fourth order equation 



{Ul - E){Ur -E)- ut^ 



1 (PUl 



72 dX2 



= 0.(2.4) 



Here 7 3> 1 is the semiclassical parameter which is determined by the ratio of the characteristic potential scale over 
the zero oscillation energy (i.e. as above 7 = miloao/^' where m is a mass of a particle, ag is a characteristic length 
of the problem, e.g. the tunneling distance, Q.q is a characteristic frequency, e.g. the oscillation frequency around the 
potential minimum). 

Luckily the equation (2.4) admits semiclassical solutions by Fedoryuk method [8] - [10] since the coefficients at the 
n-th order derivatives proportional to 7~", and therefore are small. Besides in the vicinity of the crossing point X = Q 
the diabatic potentials (2.1) can be replaced by the linear ones counted from the barrier top U"^ 

UL/R{X) = U*±fX, (2.5) 

and eventually the equation (2.4) can be presented into a more compact and simple form 

+^ -C/i2]eL=0, (2.6) 

where a = U"^ — E. 

Four roots of the characteristic polynomial of (2.4) or (2.6) 



F{\ X) = - ^\Ul + Ur- 2E)\^ - 27'^A + 7^ 



{Ul - E){Ur -E)- U?^ 



1 d^UL 



72 dX^ 



(2.7) 



determine the four fundamental solutions to (2.6) 

yj = ifX^ + Ul^)-^'^ exp (^j \j{X)dX^ , j = 1, 2, 3, 4 . (2.8) 
The solutions (2.8) can be visualized as a motion with imaginary momenta in the upper and lower adiabatic potentials 

C/± ^\{Ul + Un) ± \[{Ul - Unf + ^Ul^Y'^ . (2.9) 

As it was mentioned above in the vicinity of the crossing point one can replace (2.4) by (2.6). In the latter equation 
the coefficient at the first order derivative is small (a 7~^), and by the substitution 

= exp(Ki,2X)$i'2 , (2.10) 

where 

'^1,2 = ±7v/^ (1 ± , (2.11) 

and (5 is a first order correction (see [13]) 5 = {"^f /A^Ja). Therefore the equation (2.6) is reduced to two independent 
Weber equations with the known fundamental solutions [11] 



{e^} = |exp(±7v^X)i?_, (^) ' ^) ' exp(±7V^^)i?-i-. (^) ^ ^) } ' ^^"-^^^ 



where v = (7/722/4/Y^) is referred traditionally as the Massey parameter, and in fact it controls the main features 
of the behavior. The corrections to the indices of the parabolic cylinder functions D and to the arguments of these 
functions can be found from (2.11) and have been calculated in [13]. 

At the next step we should perform the asymptotically smooth matching of the solutions (2.8) and (2.12). The 
whole analysis can be brought into a more elegant form by introducing connection matrices which link on the complex 
plane the semiclassical solutions to the Schrodinger equation for the exact potential of the problem under study (e.g. 
(2.1) for our case) and the exact solutions of the so-called comparison equation (in our case (2.6)) which is valid near 
the crossing point. The explicit calculations of the connection matrices are rather involved since the LZ problem is 
characterized by the four fundamental solutions to the left and to the right regions with respect to turning or crossing 
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points. Therefore the connection matrices, we are looking for, arc 4x4 matrices. Although the generalization for our 
case of the known already 2x2 connection matrices (see e.g., [12], and for more recent references our publication [13]) 
is straightforward, it deserves some precaution as it implies quite different procedures for the energy, more accurately 
for E/^ smaller (the tunneling region), larger (the over-barrier region), or of the order (the intermediate region) of 
the potential barrier, i.e. U'^ — U\2- 
Indeed, in the case 



-«U* 
7 



U- 



12 , 



(2.13) 



the region near the crossing point is forbidden for the both adiabatic potentials. However, four real-valued turning 
points of the lower adiabatic potential are far enough from the crossing point. The upper adiabatic potential in 
this case is also higher than S/7, and therefore for the instanton approach there are two imaginary turning points 
which characterize the motion in the inverted upper adiabatic potential. Thus for the tunneling region we have four 
real- valued and two pure imaginary turning points. 

In the over-barrier energy region, when the energy is larger than the upper adiabatic potential minimum, i.e. 



E 
7 



(2.14) 



the whole region for the both potentials, is accessible for the classical motion. Therefore there are four real-valued 
turning points (two for the lower and two for the upper adiabatic potentials). Besides there are two imaginary 
turning points corresponding to the quantum over-barrier reflection for the lower adiabatic potential. Finally in the 
intermediate energy region, i.e. for 



U* + Ui2>->U*- U12 , 
7 



(2.15) 



there are two real-valued and four imaginary turning points. 

The tunneling path is one central point to be considered within the instanton method, and the determination of the 
tunneling trajectory (or trajectories) is, in a general case, a nontrivial task. However for our model ID potential (1.1) 
in the symmetrical case the extremal action trajectory consists from so-called kink and anti-kink pairs corresponding 
L ^ R and R ^ L transitions, and the action for every part (i.e. kink or anti-kink) is W*. More or less qualitatively 
the same is the tunneling path for a small potential asymmetry. However, when the asymmetry is larger than the 
tunneling splitting in the symmetric double- well potential, there is only one classical trajectory starting from the less 
deep well (say L) which does not reach the more deep R minimum and comes back to L. Thus in this case the pair 
kink - anti-kink forms a single so-called bounce trajectory with the action 2W*. We will explore this issue in more 
details in what follows. 

The double-well shape of the lower adiabatic potential and the influence of the upper adiabatic potential require 
that to find the solutions one has to take into account at least two instanton trajectories with the energies E = 
and E = 7^^. Following the strategy, described above, one has to match smoothly the semiclassical (e.g., instanton) 
solutions known in the remote from the crossing point {X = 0) region with the solutions of the more simple comparison 
equation which is valid in the vicinity of the crossing point. This matching should be performed asymptotically, i.e. 
at small \X\ but for large enough 

Now we are in the position to find all needed connection matrices. In the tunneling region (2.13) for every well [L or 
R) there exist increasing and decreasing exponentially real-valued solutions to the Schrodinger equation. The solutions 
are matched at the crossing point, therefore they are linked by the real- valued 4x4 connection matrix which should 
have two 2x2 blocks linking the increasing (decreasing) diabatic solution in the L-well with the dccrciasing (increasing) 
diabatic solution in the i?-well, in the agreement with the standard Landau scheme of the tunneling transitions [1]. 
Omitting a large amount of tedious algebra we can represent the connection matrix linking the "asymptotic" (i.e. 
in the left/right (L, R) wells and for the upper/lower {+, — ) adiabatic potentials) solutions in the tunneling energy 
region in the following form 



U, 



(2.16) 
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Here tJc is the 4x4 connection matrix at the crossing point, which in the tunneUng region has the following form 



P 



cos(7r!/) 





(sin^(7ri/))/p 
cos(7rz/) 






- cos(7ri^) 

P 




— cos(7ri^) 



{sin^ {ttu)) / p 



(2.17) 



where we designated 



P = 



2t t exp(-2x) 

(+) 



(2.18) 



and X = (^^/2) — (1/2) (ly — (1/2)) In z/. The matrices Mcf^' and Afi ' arc the 2x2 connection matrices at the 
corresponding linear turning points, which are determined by the phase shifts at these points 



Mi-) = 



1 -i \ 
-{i/2) (1/2) j 



(2.19) 



and Mi^) is the matrix Hermitian conjugated to (2.19). The L'^^jj^ and matrices are called shift matrices, and 
those are related to the variations of the coefficients of increasing and decaying semiclassical solutions in the regions 

between the turning points {Fc is the shift matrix when one moves from the crossing to the turning point in classically 

forbidden region, and L'^^jj^ are the shift matrices in the classically accessible regions). Explicitely we get 



P _( eM-W%/2) 
" ' exp(7T4^^/2) 



(2.20) 



Here 7 is the semiclassical parameter, and is the action in the lower adiabatic potential barrier. Finally the 
structure of the shift matrices -C'^jjj is 



f (c) 



/ exp(z7VF2/^) 

V exp(-nW^*/^) 



(2.21) 



where W^/i? action calculated by the integration between the turning points. We depicted the corresponding 

trajectories and matrices in Fig. 2 

The same manner can be treated the over-barrier region (2.14) (see Fig. 3). In this case the crossing point is in 
the classically accessible region for the both potentials. The fundamental diabatic solutions can be represented as 
the waves propagating in the opposite directions, and the complex-valued connection matrix has as it was for the 
tunneling region 2x2 block structure, where the blocks link the waves in the L and in the R wells propagating in the 
same direction. Specifically the corresponding connection matrix at the crossing point U'^ 



s exp(— i(^) — exp(— TTZ/) 

sexp(i(/)) — exp(— TTJ/) 

exp(— TTJ^) sexp(— i(/)) 

exp(— TTi/) sexp{i<p) 



(2.22) 



(where we denoted s = \/l — exp(— 27r:/), = argr(— ii^) -I- 3(2%), and x ~ ^(V2)((7r/4) -I- — hiu)) + (l/4)(7ri/ -|- 
Inj/)) should be multiplied by two blocks: the block from the left gives the contribution at the turning point and 
includes the shift matrix to the crossing point in L and in R wells of the lower adiabatic potential; the right block is 
related to the turning point and to the shift matrix to the crossing point in the upper one-well adiabatic potential. 
Thus finally in the over-barrier region we get 



L 








(2.23) 



Here we used the same notations as it was above for the tunneling region, and besides the matrices M(='=) are transposed 
with respect to the matrices Mc"^^ given in (2.19), and the new shift matrix L is 

exp(-i7W72) 

exp(i7VF*/2) 



(2.24) 
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(remind that W* is the action in the upper adiabatic potential). Combining altogether (2.23), (2.22), (2.24), and 
(2.19) one can trivially find the full connection matrix for the over-barrier energy region (2.14). 

More tricky task is to calculate the connection matrix in the intermediate energy region (2.15). In this region the 
crossing point is close to the internal linear turning points of the diabatic potentials. Therefore the two fundamental 
diabatic solutions are in the classically accessible region and two others are in the forbidden region (see Fig. 3 for 
the illustration) . Nevertheless even in this case the connection matrix has the 2x2 block structure, but these blocks 
determine the transitions between the adiabatic states (unlike the tunneling or the over-barrier regions, where the 
connection matrices (2.16), (2.23) link the diabatic states). To treat this kind of problems we have developed recently 
[13], [14] the semiclassical method for level quantization in the vicinity of the diabatic potential crossing point. The 
method enables to find all four exponentially increasing or decreasing solutions to the Schrodinger equation for an 
arbitrary shape of the crossing diabatic potentials, i.e. for any combination of the 1-st and 2-d order turning points, 
and of the crossing point. In this paper we generalize this method [13], [14] to study coherent-incoherent tunneling 
regimes in an asymmetric double well potential. 

Following the same line as above we first present the general structure of the connection matrix in the intermediate 
energy region 



(+) 







(2.25) 



The matrices 



have been introduced in [15] for the imaginary turning points (in the case tinder consideration 
the both adiabatic potentials have these turning points in the energy window ]a| < U12) and they have a form 



m: 



(+) 



1 

(*/2)exp(-7M^±) 1 



where Mj_ is the matrix Hermitian conjugated to M±^\ and are so-called Euclidean actions in the turned over 
upper and lower adiabatic potentials 



^; 91,2 =7ttvC^12±Q!, . 

7 4/ 



(2.26) 



Since the matrices M^j^-* become the unitary ones at a > [/12 and at a < — 1/12, the intermediate region connection 
matrix (2.25) matches continuously the the connection matrices (2.16) and (2.23) in the tunneling and over-barrier 

regions, respectively. 

The connection matrix in the intermediate energy region can be calculated using the known Weber function asymp- 
totic expansions for large complex indices [16], [17]. Combining together these asymptotic expansions and all matrices 
entering (2.25) defined already above, we find at the crossing point, the matrix TJ" is 



U' = 



(V27r/r(g*))exp(-2x(9*)) 

(r(9)/^)exp(2x(g))(l - exp(-27rg2)cos2(7r9i)) 
exp ( — 27r(72 ) cos (tt^i ) 

exp(— 27rg'2) cos(7rgi) 



(2.27) 



— exp(— 27792) cos(7rgi) 

— exp(— 27792) cos(7rgi) 
(x/2^/r(g))exp(2x(9)) 

(r(9*)/y2^) exp(2x(9*))(l - exp(-27r92) cos2(7r9i)) 



where as above we introduced the following abridged notations 



7Mi2\/^*12 ± OL * 

9 = Qi + 192 ; 91,2 = ; 9 = 91 - *92 , 



(2.28) 



and, besides. 



X = xi + «X2 ; 2x1 = 91 - ( 91 - 2 ) In l^l + mi ■ 



(2.29) 
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and analogously 



The phase factor is defined as 



2X2 = g2-«2ln|(z|-<^(gi-^ ) . (2.30) 



Now the full connection matrix in the intermediate energy region can be found easily simply collecting the given above 
expressions. 

III. RESULTS AND DISCUSSION 

Our purpose in this section is to study how found in [2] coherent - incoherent tunneling relationships shortly 
described in the introduction section I for an isolated double- well potential, and in particularly the criterion (1.4) and 
the dephasing time T2 (1.3) should be modified for more realistic situations when there is a finite adiabatic coupling 
between the diabatic potentials forming the asymmetric double-well lower adiabatic potential, and the one-well upper 
adiabatic potential. 

However to investigate this problem first we should derive the quantization rules for the crossing diabatic potentials. 
It can be done using presented above the connection matrices. In spite of the fact that instanton trajectories are 
rather simple objects, and can be relatively easy found analytically, calculations of the quantization rules within the 
instanton approach are rather intricate and require the knowledge of all connection matrices, we have calculated in 
the previous section. To apply this machinery within the instanton approach, the quantization rule can be formulated 
as a condition that the amplitudes of exponentially increasing solutions must be vanished. In terms of the matrix 
elements of the connection matrix this condition is 

TO22m33 - m23m32 = . (3.1) 

If made above assumptions are granted one can easily write down the Bohr - Sommerfeld [1] quantization equations 
applying shown in the Figs. 2, 3 the connection- and shift matrices (for the details of the calculation method, see 

[13], [14]), and we end up with 

tan(7W£) tan(7W^^) = 1 exp(27W^,*) , (3.2) 

where is the action in the classically forbidden region in between the turning points, and W^jj^ are the coordinate 
independent actions inside of the L (respectively R) well. This equation (3.2) can be solved to find energy levels in 
the wells. 

Applying the same procedure to the over-barrier region (2.14) we find from (2.26) 

(1 - exp(-27ri^)) cos(7(VK£ + W},) - 0) cos(7T4^* +4))+ (3.3) 
exp(-27r!/) cos (wi + cos (w}^ + ^ 

In the diabatic limit [u — > 0) one get from (3.3) 

cos (^7 (wi + ^) ) cos (^^(w*^ + ^y^=Q, (3.4) 

and therefore two independent quantization conditions 

On the other hand in the adiabatic limit, i.e. at ^ 00 and ^ ^ we have 

cos(7(W^£ + W^)) cos(7VF*) = , (3.6) 
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and therefore 



^{Wl + W^) = n(n+^) ; iW*)=^(no + l] ■ (3.7) 



And to conclude this part and to span a wide range of possibihties, the quantization condition in the intermediate 
energy region derived from the corresponding connection matrix can be represented in the following form 

cos(7(W^£ + W^) -</.)= ,;'^P^~7^^ ^ cMjiWl W^)) . (3.8) 

Vl + exp(-7rg2) 

Now it seems appropriate to take a fresh look at the results presented above. What can we learn from the performed 
calculations? First we can go one step further to analyze the phase factors calculated above. In our system (two 
crossing diabatic potentials) there are two types of phases. The first phase factor occurs, since the tunneling results 
in the phase shift related to the change of eigenvalues. In own turn it leads to a certain kind of one- well phase 
(T2) relaxation. The physical argument leading to T2 relaxation at the tunneling in the asymmetrical double well 
potentials may be rationalized as follows. The fact is that reflected from the barrier waves acquire non trivial phase 
factor. The phenomenom is related to interference of incident, reflected and transmitted waves. One can look to this 
phase factor from a slightly different point of view, since tunneling results in the phase shift related to the change of 
eigenvalues. The quantization rules can be rewritten in the form which includes some integer numbers numerating 
and an exponentially small phase shift due to the existence of the barrier between two wells. 

The second phase shift occurs in our case due to non-adiabatic behavior. Indeed the LZ case (even for the same 
asymmetric double well shape of the lower adiabatic potential) is quite different not only quantitatively due to coupling 
with the upper adiabatic potential but also qualitatively, since a novel and fundamental quantum effect will occur. 
Namely, in addition to the described above tunneling phase (existing even in an isolated double well potential) a 
quantum mechanical wave function acquires upon a cyclic evolution some geometrical or Berry phase factor [21] - 
[22]. What is most characteristic for the concept of Berry phase is the existence of a continuous parameter space in 
which the state of the system can travel on a closed path. In our case the phase is determined by the non-adiabatic 
interaction. Coherent or incoherent kind of behavior for crossing diabatic potentials crucially depends on a quite tricky 
interplay between the both (i.e. tunneling and Berry) phase factors. Two new results which have emanated from our 
study of these phenomena, is our main motivation for presenting this paper. The semiclassical wave functions of the 
bound states are linear combinations of (2.16), or (2.23), or (2.25) (respectively for the tunneling, over-barrier, and 
intermediate energy regions), which can be determined, provided we know the eigenvalues. The quantization rules 
(namely, conditions that coefficients at exponentially increasing and ingoing from the both infinities jXj ±00 waves 
are zero), and the wave function normalizations, define uniquely these linear combinations. The shift matrices have 
always the same form as (2.21), where one has to insert the action between a given point X and the nearest turning 
(or crossing) point Xk- For the upper and the lower adiabatic potentials the action reads as 



W±{Xk,X) = £ dxJ2 



(3.9) 



We have illustrated the general connection matrix scheme in the Figs. 2 and 3. 

Let us consider a general example describing two non-symmetric potentials crossing at X = point (2.1). When 
the parameter h entering the potential (2.1) is varied from 1 to 00, we recover the two known in the literature limiting 
cases, and come from two identical parabolic potentials to the case one-well and linear diabatic potentials crossing. 
This kind of the diabatic potentials crossing leads to the lower adiabatic potential in the form investigated in [2] , and 
has qualitatively the same features as the model potential (1.1). If one neglects for a moment the upper adiabatic 
potential, aiming to study crossover behavior from coherent to incoherent tunneling upon increase of the parameter 
6, the larger is this parameter 6, the larger will be the density of final states. The criterion for coherent-incoherent 
crossover behavior found in [2] based on comparison of the transition matrix elements and the inter level spacings 
in the final state. The analogous criterion should hold for LZ level crossing problem, however in the latter case the 
tunneling transition matrix elements has to be multiplied by the adiabatic factor. Therefore the coherent - incoherent 
tunneling crossover region moves towards the more dense density of final states, and the larger Ui2 is the smaller 
will be the region for incoherent tunneling. For the sake of the skeptical reader it is worth to emphasize that the 
tunneling matrix element dependence on the Massey parameter v found above is valid for an arbitrary magnitude of 
the adiabatic coupling (cf. the recent publication [38] where this matrix element has been calculated in the frame 
work of the perturbation theory, and only in the adiabatic and diabatic limits). 

Owing to the non-adiabatic behavior of the system the tunneling matrix element Hlr is renormalized by the 
adiabatic factor. In the tunneling region from (3.2) we find this renormalization as 

Hlr ^ Hlrv{v) , (3.10) 
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where the function p{u) (2.18) is associated with the transition amplitudes between the diabatic potentials in the 

crossing region. 

This renormalization tunneling factor varies from to 1 upon increasing of the Massey parameter u. Again as we 
have found for the isolated double- well potential, in the limit 

Hlrp{u)^Ar, (3.11) 

the spectrum consists of the set of the tunneling doublets and L — R transitions are coherent ones. This criterion 
(3.11) replaces (1.4) for our case of the finite adiabatic coupling U12. 

Quite different situation occurs for the excited states. In the diabatic limit, the transition matrix element is increased 
with the Massey parameter v, and therefore at a given b value, the system moves to more incoherent behavior. In 
the adiabatic limit, the transition matrix element is exponentially small, and coherence of the inter-well transitions 
should be restored. However, since the matrix elements are oscillating functions of C/12 for the intermediate range of 
this coiipling (t/12) coherent - incoherent tunneling rates are also non-monotonically varying functions. To illustrate 
it let us study dynamics of the initial state of the system, prepared somehow in the ground vibrational state of the 
upper adiabatic potential 



f/2 + X2 I 1 + "—^X 



'12 



6-1 



(3.12) 



Evidently the wave function ijA^'' of this state should be close to the harmonic oscillator function n = 0, with its 
minimum at X = 0, and its eigen-state Eq dependent of U12. The real part of the i?o gives the oscillator vibration 
frequency wqj and the imaginary part determines the decay rate To of the state. 

The spectral expansion of ijj^^^ over the diabatic state eigen-functions {V'n} can be found by the method proposed 
long ago by Zeldovich [27] (see also [2], where the method has been applied to find semiclassical solutions of the 
Schrodinger equation for ID asymmetric double well potential) well adapted for quasi-stationary state wave function 
expansion over continuum spectrum functions. Note that for 6 ^ 1 all the states of the upper adiabatic potential 
(3.12) are placed in the same energy range that the ii-well excited states. Since, as it is well known [1], the harmonic 
oscillator wave functions for n 3> 1 coincide with the semiclassical ones, we can represent the i?-well wave functions 
in the vicinity of Eq in the following form 

f ^(^.W4"'(A-).|.Y|.A„ 

where Aq = h/muo is the de Broglie wave length, corresponding the ground state oscillator wave function, ipn^ are 
harmonic oscillator eigen - functions, A;„ = ^J^mEnlTi is the wave vector, and the phase factor v'(fcn) is defined as 



The entering (3.13) amplitudes A(fc„) are determined by the condition that the probability density flow from the quasi- 
stationary state to infinity should be constant (in fact this condition plays the role of the normalization condition for 
the quasi-stationary states [27]): 



We have shown in the Fig. 4 (for a fixed value of the potential shape controlling parameter 6 l^- 1 and various coupling 
strengths Uyi) the spectral density expansion 

S{E)=Y,\<i,\i,^>\H{E-En). (3.16) 

n 

The spectrum of final states for 6 ^ 1 is a discrete one (although dense), and the envelope of the spectrum has a 
Lorentzian shape with a width Fq, which is determined by the non-adiabatic transition matrix element. The latter 
quantity has an oscillating dependence on C/12. Since the final state spectrum at 6 » 1 has only weak dependence on the 
non-adiabatic coupling l]\2 we have almost constant spectral distribution but the number of the final states (relevant 
for the transition) is determined by Fq, and it oscillates with Uyi- Analogously to the adiabatic transitions considered 
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in [2], one can formulate the non- reversibility criterion for non-adiabatic transitions. Indeed coherent oscillations turn 
into exponential decay when sufficiently large number of the final states occurring under the Lorentzian envelope, i.e., 
when 



It is worth noting that unlike the adiabatic transitions [2] , when Fq is determined uniquely by the amplitudes of the 
inter - well transitions, in our case for the quasi-stationary states of the upper adiabatic potential, Fq depends on the 
rate of the non-adiabatic transitions. To calculate the rate one has to solve first the eigenvalue problem. 

The spectral distribution (3.16) fully characterizes the state time evolution. In particularly relevant quantity 
is so-called survival probability P{t) reads as 



and it is presented in the Fig. 5. We see that there are a number of the coupling strengths U12 regions with the 
non-monotonous behavior. The coherent dynamics holds for the energy regions not too close to the L well levels, when 
only few R states coupled to the state i'it). In the case the energy of the state is in the resonance with a certain level 
of the L well, the transition matrix elements arc enhanced. In this case the ip{t) evolution resembles decays of the 
quasi-stationary states (see e.g., [13], [14]). It is convenient and instructive to illustrate the non-monotonous behavior 
by averaging of the surviving probability over the period of the initial state recovering. Evidently this average value 
< P >T equals 1/2 for the coherent oscillations. The dependence of < P >t on U12 is presented in the Fig. 6. 
Note that in the both (adiabatic and diabatic) limits we have coherent oscillations, as it should be (although over 
very different time scales), while in the intermediate energy region the coherent evolution is multiply destroyed in the 
vicinity of the resonances. 



To conclude let us comment first on our motivation. In principle potentials with two stable equilibrium config- 
urations are widely used in chemistry and physics, to describe molecular spectroscopy data. Analyzing these data 
one must distinguish two typos of states which require a set of different theoretical and experimental methods, each 
one with specific strengths and weaknesses on certain scales of length and time. Experimentally low-energy states 
localized near the minima of such potentials, are studied by vibrational spectroscopy methods. These low-energy 
states can be characterized by well defined quantum numbers describing the normal vibration excitations. Quite 
different approaches are necessary to use to study highly excited states situated near the potential barrier top. Just 
these states determine the probability of thermo-activated molecular transitions. These phenomena are intrinsically 
statistical ones, and not surprisingly theoretical descriptions of these states assume usually their ergodicity. This 
ergodic behavior can be easily understood since the excited states near the barrier top have so high density that even 
very small coupling to an environment (thermal reservoirs) can provide say fast mixing and thermalization of the 
states. 

However applying these approaches (and the model potentials) to real chemical dynamic problems of low- 
temperature reactions and transitions of relatively small molecules or atomic clusters (attracting much attention in 
relations with chemical reactions in upper Earth atmosphere layers, and high precision laser spectroscopy techniques), 
one should take care whether these two relevant regions of energy are not overlapped. Measurements of molecules 
with two stable configurations performed in the temperature interval (10 — 20 K) low enough to provide that for the 
measurement time dephasing or relaxation processes arc not essential [28] , [29] , [30] , demonstrated tunneling doublets 
dependence on well defined vibrational excitations. Thus this low-temperature behavior can be attributed with the 
coherent tunneling, and the advent of ultrafast lasers has provided physical chemists with a tool for studying these 
systems under nonequlibrium conditions. 

On the other hand there are also numerous examples (see e.g. [31], [32], [33], [34]) clearly showing exponential 
(incoherent) decay of an initially prepared seemingly equilibrium configuration. Thus experimental data signal that 
two distinct dynamic regimes exist in bistable molecular systems. Moreover, barrier heights and potential well 
asymmetries in the systems manifesting different dynamic behaviors, are quite similar by their magnitudes (of the 
same order). A question of primary importance is the understanding of how these two tunneling regimes (coherent 
or incoherent) depend on more subtle specific features of the potential energy profiles than merely energies of the 
characteristic points (like for example, our model potential (1.1), (2.1)). One qualitative answer to this question has 
been done long ago by Jortncr and Bixon [35]. In this paper the authors have formulated the irreversibility criterion. 
According to the criterion, coherent tunneling should be destroyed when the density of final states is so high that 



Fo > AEn . 



(3.17) 




2 



(3.18) 



IV. CONCLUSION. 
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typical inter-level spacings become smaller than characteristic transition matrix elements. Our aim in this paper is 
to formulate quantitatively the analogous criterion. One methodical comment seems in order here. Usual technique 
to analyze radiationless transitions is based on the perturbation theory in the adiabatic representation, and the non- 
adiabatic coupling operator is treated as a perturbation (see, e.g., [41], [42]). It is easy to understand, however, that 
this kind of the adiabatic perturbation theory is equivalent to isolated two level systems approach, which is valid only 
when the level displacements (due to tunneling or non-adiabatic transitions) are smaller than the inter-level spacing. 
Analogously for quasi-stationary states the adiabatic perturbation theory works when the level broadening (or decay 
rate) is smaller that the level spacing. Clearly it is not the case for the intermediate energy region we have studied 
in our paper. Note also that our connection matrix approach can be also formulated in terms of the Liouville - von 
Neumann equation for the density matrix, where so-called relaxational matrix should be chosen phenomenologically 
to mimic decay rate dependences on the energy and on the Massey parameter. But anyway to find those one has to 
solve the Schrodinger equation for the potential under consideration. 

Our model appears to be a simplest one demonstrating that relatively small variation of the adiabatic coupling 
(at a level spacing scale which is small in comparison with potential barrier heights or well asymmetries) enables 
to change qualitatively dynamic behavior. Therefore, we conclude that dynamic irreversibility in the systems under 
investigations crucially depends on the final states density (and not on potential energy profiles directly). We believe 
we are the first to explicitely address this issue. To illustrate these phenomena we investigated coherent and incoherent 
tunneling in the conditions of crossing diabatic potentials. As a result of the coupling and level crossing avoiding, we 
get the asymmetric double-well lower adiabatic potential with a variable shape depending on a value of a continuous 
parameter b (which describes in the limit 6=1 two identical parabolic diabatic potential crossing and in the limit 
b ^ 00 one- well and linear diabatic potentials crossing). The doublet structure of levels (generic for double- well 
potentials) is remained valid as long as the renormalized by the adiabatic coupling transition matrix element Hlr 
(or tunneling splitting) is smaller than characteristic inter-level spacings Ajj We calculated the non-adiabatic factor, 
and found in the diabatic limit (U12 0) H^n goes to zero, and in the adiabatic limit {U12 — > 00) the tmmeling 
transitions do not depend on the upper potential. In the over-barrier energy region Hlr is an oscillating function 
of U12, due to the resonances between the states in the lower and in the upper adiabatic potentials. In the case 
Hlr > A.R, any level from the shallow L-well is coupled by the tunneling to several levels in the i?-well, and the 
transitions lose their coherence. 

In an apparently unrelated development researches studying the problem of intermolecular energy redistribution 
discovered purely quantum energy flow between modes which would be otherwise uncoupled. The mechanism for such 
classically forbidden energy flow between degenerate vibrational modes arose from non-adiabatic couplings involving 
a sequence of intermediate states. In the various existing in the nature molecular systems the non-adiabatic coupling 
strengths U12 can have fairly diff'erent magnitudes, thus we anticipate realizations of the both kinds of dynamic 
behavior, coherent and incoherent ones. Our model of the non-adiabatic transitions from the initially prepared quasi- 
stationary state ^j{t) can be directly confronted to experimental data on super fast non-linear optical spectroscopy (see 
e.g., the monographes [39], [40]). The developed in this area technique allows to prepare a given initial quasi-stationary 
state by a suitable optical pumping pulse shape. 
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Figure Captions. 

Fig. 1 

Crossing asymmetric parabolic diabatic potentials (adiabatic potentials for U12 = 0.5) are shown by the dashed 
lines): 

(a) bound initial and final states (6 = 3); 

(b) bound initial and decay final states {b = 00). 
Fig. 2 

Connection matrices for the tunneling energy region: 

(a) in the WKB approach, where the trajectory has 4 linear turning points, and one crossing point. mJ^^ are 
the connection matrices for the isolated linear turning points, L^'j^jj^ are the shift matrices (2.21) in the classically 
accessible regions, is the shift matrix (2.20, and Uc is the connection matrix at the crossing point; 

(b) in the instanton method, where the trajectory E = passes through the second order turning point (the L 

well potential minimum). In the L well the WKB connection matrix should be replaced by the energy dependent 

(2) 

connection matrix Mj^ at the second order turning point [15]), c;orrcsponding to kink-anti-kink pair. The same 

(2) 

manner the connection matrix for the R well corresponds to so-called bounce [25] . 
Fig. 3 

Connection matrices for the over-barrier and intermediate energy regions (in the latter region the connection 
matrices for the imaginary turning points are not shown). 
Fig. 4 

Spectral distributions for the initial state wave function (chosen as the ground state of the upper adiabatic potential) 
over the eigenfunctions of the model potential (1.1): b = 1500, 7 = 12, and U12 = 0.09 , 0.15, 0.21 , 0.28 , 0.40 for (a) 
- (e) figures respectively. 

Fig. 5 

Survival probability for the same as in Fig. 4 initial state. 

(a) - solid line corresponds to the Fig. 4a; dashed line - to the Fig. 4b; 

(b) solid line corresponds to the Fig. 4c, dashed line - to the Fig. 4d. 
Fig. 6 

Averaged over the recovering period survival probability shown in Fig. 5. 
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